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Abstract 



We study the cosmology of the Brans-Dicke theory with perfect fluid type 
matter. In our previous work, we found exact solutions for any Brans-Dicke 
parameter uj and for general parameter 7 of equation of state. In this paper we 
further study the cosmology of these solutions by analyzing them according 
to their asymptotic behaviors. The cosmology is classified into 19 phases 
according to the values of 7 and ui. The effect of the cosmological constant 
to the Brans-Dicke theory is a particular case of our model. We give plot of 
time evolution of the scale factor by numerical investigations. We also give a 
comparison of the solutions for the theories with and without matter. 
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I. INTRODUCTION 



Recent developments of the string theory suggest that in a regime of Planck length 
curvature, quantum fluctuation is very large so that string coupling becomes large and 
consequently the fundamental string degrees of freedom are not weakly coupled good ones 
IJ. Instead, solitonic degrees of freedom like p-brane or D-p-brane [JJ are more important. 
Therefore it is a very interesting question to ask what is the effect of these new degrees of 
freedom to the space time structure especially whether including these degrees of freedom 
resolve the initial singularity, which is a problem in standard general relativity. The new 
gravity theory that can deal with such new degree of freedom should be a deformation of 
standard general relativity so that in a certain limit it should be reduced to the standard 
Einstein theory. The Brans-Dicke theory M is a generic deformation of the general relativity 
allowing variable gravity coupling. Therefore whatever is the motivation to modify the 
Einstein theory, the Brans-Dicke theory is the first one to be considered. As an example, 
low energy limit of the string theory contains the Brans-Dicke theory with a fine tuned 
deformation parameter (uj=-1). 

Without knowing the exact theory of the p-brane cosmology, the best guess is that it 
should be a Brans-Dicke theory with matters. In fact there is a some evidence for this ||, 
where it is found that the natural metric that couples to the p-brane is the Einstein metric 
multiplied by certain power of dilaton field. In terms of this new metric, the action that 
gives the p-brane solution becomes Brans-Dicke action with definite deformation parameter 
u depending on p. Using this action, Rama |J recently argued that the gas of solitonic 
p-brane || treated as a perfect fluid type matter in a Brans-Dicke theory can resolve the 
initial singularity without any explicit solution. In a previous papers @,@]> we have studied 
this model and found exact cosmological solutions for any Brans-Dicke parameter uj and for 
general equation of state and classify the cosmology of the solutions according to the range 
of parameters involved. 

In this paper we further study the cosmology of these solutions by analyzing them ac- 
cording to their asymptotic behaviors. The cosmology is classified into 19 phases according 
to the values of 7 and uj. The effect of the cosmological constant to the Brans-Dicke theory 
is a particular case of our model. We give plot of time evolution of the scale factor by 
numerical investigations. We also give a comparison of the solutions for the theories with 
and without matter. 

The rest of this paper is organized as follows. In section II, we set up the notation and 
review the result of our previous results H^fl- In section III, we describe two new phases 
which were not mentioned in |7j] and classify the cosmology into 19 phases. By the numerical 
work as well as the analytical method, the behaviors of the scale factor is presented explicitly 
by figures. In section IV, we summarize and conclude with some discussions. 

II. BRANS-DICKE COSMOLOGY WITH MATTER 

First, we briefly review our earlier work 0. We consider the Brans-Dicke theory and 
analyze the evolution of the D dimensional homogeneous isotropic universe with the perfect 
fluid type matter. The action is given by 
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S = J d D xy/=ge~* [K - ud^(j)\ + S n 



(1) 



where <fi is the dilaton field and S m is the matter part of the action. Here we assume that 
the matter has no dilaton coupling. 

Let's assume that the matter can be treated as a perfect fluid with the equation of state 



p = 7 p,7 < 1. 

Therefore our starting point is the equation of the Brans-Dicke theory P,|TU 



(2) 



K 



g -^n = ^ v + uj{d^d^ 



(<90) 2 } 



+{-d (t d v <j> + d^d v <\) + g^V 2 ^ - g^(d(fr) 2 }, 
® = n-2ujV 2 <p + uj(d<p) 2 



(3) 



where <fi is the dilaton and T> means a covariant derivative. 1Z is the curvature scalar and 
the metric is given as the follows: 



ds 2 = -jfdt 2 + e 2a{t) 5 ij dx i dx j (i, j = 1, 2, ■ ■ ■ , D - 1), 



(4) 



where e a ^(= a(t)) is the scale factor and J\f is the (constant) lapse function. 
The energy-momentum tensor is given by 

Tp> = pg^u + (p + p)U^U u 



(5) 



where is the fluid velocity. The hydrostatic equilibrium condition of energy-momentum 
conservation is 

p+(D-l)(p + p)a = 0. (6) 
Using the equation of state, eq.(0), with a free parameter 7, we get the solution 

p = poe-^-W+T)", (7) 

where po is a real number. Our goal is to study how the metric variables change their 
behavior for various values of 7 and to. Now, since we consider only the time dependence, 
the action can be brought to the following form 



S 



dt e^- 1 ^ 



1 



{ — (D — 2)(D — l)a 2 + 2(D - l)d0 + cj0 2 } 



-VMp e 



-(D-l)(l+j)a+cf> 



(8) 



where we eliminated p and p by eq.(|6|). After getting the constraint equation by varying 
over the constant lapse function, A/", we can set it to be 1. 
Now, introducing a new time variable r by 



dt 



,(D-l)a- 



(9) 
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and the new variables 



X = _-[p-l)(l_ 7 )a- 
Y = a + —X, 

K 



the action can be written as 



S = dr 



{(D - 1)kY 2 + pX 2 } - VATpoe 



-2X 



(10) 



(11) 



where 



K 
V 

fl 



(D-l)(l--y) z (u-u R ), 
2{l- 1 ){u-u v )i 
A(D - 2) 



K 



0J K 



D - 2D~i + 27 

'(£-l)(l-7) 2 ' 
1 

1^7' 
D - 1 

D - 2' 



The constraint equation is given by 

— (D — 1)kY 2 + pX 2 + p e~ 2X . 



(12) 



(13) 



Note that for D > 2 and 7 < 1, the sign of k is determined by that of uj — uj k and the sign 
of p is determined by that of uj — oj-\ and k. 
The equations of motion are simply 








Y, 

X - ^e- 2X . 



(14) 



When p = 0, the situation is that of the string cosmology discussed first in and the 
solution for X is X = cr. One can easily show that this solution has two disconnected 
branches in terms of the original time £; one is inflation-type and the other is FRW-type. 
If po 7^ 0, the asymptotic behavior of X is X ~ c | r | as we will see later. In other word, 
the behavior of cosmology at p 7^ is not continuously connected to that of cosmology at 
po — in r — ► —00 region. Some of the new aspects of the cosmology due to the presence 
of the matter come from this discontinuity. If po is a negative constant, then the solution 
oscillates in time, leading to unphysical solution. This includes the situation where there is 
a negative cosmological constant in the Brans-Dicke theory. In this paper, therefore, we will 
consider only positive po- 

If uj is less than w_ l3 the kinetic term of the dilaton has a negative energy in Einstein 
frame. So we will consider the case where u is larger than uj_i. According to the sign of k, 
the types of solutions are very different. When k is negative, the exact solution is 



4 



X — In [- cosh(cr)], 



Y = At + B, (15) 
iftrary real i 

can determine A in terms of other parameters, 



where c, A, B and g = Jjj^ are arbitrary real constants. Using the constraint equation, we 



c . , / (D — 1)k \ k\ , . 

A=~, with 6 = \ -- — = — ; ' ' 16 

5 V A* 2 v /l+^§Ef 



If k is zero, it turns out that the solution of the equations of motion does not satisfy the 
constraint equation. If k is positive, the solution is 

q 

X = In [- | sinh(cr) | ], 

Y = -r + B. (17) 
o 



III. PHASES OF THE BRANS-DICKE THEORY 

In [0, the behaviors of the scale factor were classified by 16 phases according to u and 
7, see figure 1, and we showed that the asymptotic behaviors of scale factor a(r) and time 
£(t) at r — ► ±oo are as follows: 







T± = 


2c 




a(r) « 




F ± = 


2c 





p-lWl+w 



D-2 
L> - 1 



T sign(K){fi: + (D - 1) 7 (1 + w(l " 7))} 



_D _ 2 

'1 + w— — t =F sign(«){l + w(l - 7)} 



D - 1 

Note that the range of t is determined by the sign of T± 



(18) 



[-00, 


00) 


if 


T_ 


< < T+, 




(-00 


tf) 


if 


T_ 


< and T_|_ 


< 


(u, 


00) 


if 


T_ 


> and T + 


> 


(ti 


tf) 


if 


T + 


< < T_. 





For k > 0, i(r) and a(r) behave as 



-signer J 



q -r, e (D-l)jB 1 

(77 - 1) I r |"- r 

2(1- 7 )(l;-^) 



(19) 
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as t goes to zero, where rj = 2 + \ D ~ l n v . t( r ) i s singular at r — > if r] > 1. So for k > and 
rj > 1, the scale factor a(£) has two branches. The asymptotic form of a(t) as a function of 
t is given by 

a(t) « [T_(t - ti)] H ~ /T - at -oo, 

a(t) « [T+(t - t f )} H+/T+ at r -> oo, (20) 

where t{ (tf) is a starting (ending) point at a finite time. Eq.(|18D contains the cases where t 
starts from — oo and/or ends at oo by setting = and/or tf = 0. According to the sign of 
T± and H±/T± and the singularity at r = 0, we classified the behavior of the Brans-Dicke 
theory @. Here we summarize the result by table 1. 



phase 


sign of k 


sign of T_ 


sign of T_|_ 


range of t 


sign of i?_/T_ 


sign of 


j 








[Hi t/J 


r 




II 








{-oo, tf] 
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III 


+ 


+ 




[U,t f ] 


+ 




III+ 


+ 
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[U, oo) 
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(— oo, oo) 


+ 


+ 


V 




+ 


+ 


[U, oo) 




+ 


VI 




+ 


+ 


[U, oo) 


+ 


+ 


vir 


+ 


+ 




[U, oo) 


+ 




VII + 


+ 




+ 


(— oo, oo) 




+ 


viii- 


+ 


+ 




[ti,oo) 


+ 




VIII + 


+ 






(-oo,t/] 






ix- 


+ 


+ 




[ti, oo) 


+ 




IX + 


+ 




+ 


(— oo, oo) 




+ 


x- 


+ 


+ 




[ti, oo) 


+ 




x + 


+ 






(-OO,*/] 






xi- 


+ 


+ 




[U, oo) 


+ 




XI+ 


+ 






(-oo, tf] 




+ 



Table 1. The sign of T± determines the range of time t as following; t(r) maps the 
real line of r to (1) if T_ < < T+, (-00,00) (2) if T_ < and T+ < 0, (-00, t f ) (3) if 
T_ > and T + > 0, [t h 00) (4) if T+ < < T_, [U,t f ]. The sign of H±/T± determines the 
asymptotic behavior of scale factor a(t). 

Now, notice that not only the sign of H±/T± but also that of H±/T± — 1 is important 
because the the universe will accelerate if H±/T± — 1 > and decelerate if H±/T± — 1 < 
when t — > ±00. Therefore, we further classify the phases of cosmology accordingly. 

A. Case uj < lo k 
1. H-/T- > 1 
• For T_ > 0, the condition H_/T_ > 1 is reduced to 
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L D-2 (D- 1)7-1 

> 7^ -i 1 r. (21) 

V D-l (1- 7 )(D-1) 1 } 

If 7 < 1/(D — 1), the condition is automatically satisfied. If 7 > 1/(D — 1), inequality 
( PIP turns out to be reduced to uj > uj k , which is surprising. This means that there is 
no solution. Therefore among the regions I and VI which have X_ > 0, H > and 
uj < uj k , only / satisfies H-/T- > 1. 

If XI < 0, the condition iX_/X_ > 1 is reduced to 



D-2 (D- 1)7-1 

t/l + u; < 7^ rr 1 r, (22) 

V D-l (1- T )(X»-1)' 1 J 

whose solution is 7 > 1/(D — 1) and u; < uj k . Only /V satisfies conditions, to < uj k , 
X_ < 0, iX_ < and JX_/X_ > 1. 



2. #+/T+ > 1 
If X + > 0, the condition H + /T + > 1 implies 



,, B "2 (J -1)7-1 ,„> 

V 1 + w d3T < -(i- 7 )(d-i)' (23) 

whose solution is given by 7 < 1/(1? — 1) and < There is no region satisfying 
X+ > 0, 7 < l/(£> - 1) and uj < uj k . 

If X+ < 0, H+/T+ > 1 is reduced to 



/-, D-2 (D- 1)7-1 

whose solution is 7 > 1/(X) — 1) or w > uj k for 7 < 1/(1? — 1). There is no region 
satisfying the conditions, uj < uj k , T + < and H + /T + > 1. 



B. Case uj > uj k 
1. H-/T- > 1 
For X_ > 0, the condition 7X_/X_ > 1 is reduced to 



D-2 (D- 1)7-1 

whose solution is 7 < 1/(D — 1) and uj < uj k . Therefore there is no solution satisfying 
the conditions, uj < uj k and JX_/X_ > 1. 
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• If T_ < 0, the condition is given by 



/ D-2 (D- 1)7-1 

The solution is 7 > 1/(D — 1) or uj > uj k for 7 < 1/(D — 1). Therefore, the solution is 
summarized by uj > uj k . But in the case uj > uj k , there is no region satisfying T_ < 0. 



2. H+/T+ > 1 
For T + > 0, the condition H + /T + > 1 is reduced to 



D-2 (D- 1)7-1 

The solution is 7 < l/(D — 1) or w > uj k for 7 > 1/(D — 1). So the solution is 
summarized by w > w K like the last case. The regions III + , VII + , IX + have the 
solution satisfying uj > uj k , T + > and H + /T + > 1. 

If T + < 0, the above condition is reduced to 



D-2 (D- 1)7-1 

whose solution is 7 > 1/(D — 1) and c<j < Therefore, there is no solution satisfying 
uj < uj k , because the solution uj < uj k is inconsistent with the assumption uj > uj k . 

3. The power behavior of scale factor at r — > 
At r — > 0, a(t) is given by 

a(t) Ex I t | r (29) 

where T is given by 

r 2(1 -7)^-1^) 
(t7-1)k 

and a constant £ becomes 

2(l- 7 )(^-u)„) pT, 2(D-l) 7 (l- 7 )(^-^) 1 

£'=[5(77-1)] e L J . 

For uj > uj k and 77 > 1, the condition that T is positive is satisfied in the region uj > uj v . The 
condition T > 1 is reduced to 

(1 - 7) + \(D - 1)7 + (D - 3)}uj + (D - 2) < 0. (30) 
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This gives the following solution 

D - 2 , D ~ 3 

- < - (1 --,)( (J - lh + J - 3 ) r "> > < 31 > 
D-2 r D-3 . , 

" > - ( 1 --,)( (J - lh + J - 3 ) for "I < ~—r < 32 > 

Note that for 7 > — ;§zf , ^ K is always greater than — ^z^uT^z^ffDZgj ■ So in that case, 
there is no solution. As a result, the solution to T > 1 is given by eq.(^). This divide the 
region XI of figure 1 into two region: for < T < 1 we call this as region XI and for T > 1 
we call this as region XII. See figure 2. 

Now we summarize all possible phases in Table 2. 



phase 


sign of 

K 


sign of 


sign of 
T+ 


range of 
t 


H_/T_ 


r 

(r->0) 


H+/T+ 


/ 


- 


+ 


- 


[ti,t f ] 


H_/T_ > 1 




< H+/T+ < 1 


II 


- 


- 


- 


(-00, t f ] 


H_/T_ < 




< H + /T + < 1 


III 


+ 


+ 




[U,t f ] 


< #_/T_ < 1 






III+ 


+ 




+ 


[ti, 00) 






H + /T + > 1 


IV 






+ 


(—00, 00) 


iJ_/T_ > 1 




< H+/T+ < 1 


V 




+ 


+ 


[U, 00) 


H_/T_ < 




< H+/T+ < 1 


VI 




+ 


+ 


[U, 00) 


< H_/T_ < 1 




< H+/T+ < 1 


vii- 


+ 


+ 




[ti, 00) 


< H_/T_ < 1 


r < 




VII + 


+ 




+ 


(—00, 00) 




r < 


H+/T + > 1 


vnr 


+ 


+ 




[ti, 00) 


< H_/T^ < 1 


r < 




VIII+ 


+ 






(-00, tf] 




r < 


< 


ix- 


+ 


+ 




[ti, 00) 


< F_/T_ < 1 


< r < 1 




IX+ 


+ 




+ 


(—00, 00) 




< r < 1 


H+/T+ > 1 


x- 


+ 


+ 




[U, 00) 


< #_/T_ < 1 


< r < 1 




X+ 


+ 






(-OO,*/] 




< r < 1 


iJ + /T + < 


xi- 


+ 


+ 




[U, 00) 


< i?_/T_ < 1 


< r < 1 




XI+ 


+ 






(-00, t f ] 




< r < 1 


< H+/T+ < 1 


xii- 


+ 


+ 




[U, 00) 


< #_/T_ < 1 


r > 1 




XII+ 


+ 






(-00, t f ] 




r > 1 


< H+/T+ < 1 



Table 2. All possible phases are classified. Here phases XII and XII + are new phases. 



By the numerical work, the explicit scale factor behaviors of all phases are shown in 
following figures. Here we set D = 4. 

All phases in figure 3 have no singularity at r = and the phases II and IV have no ini- 
tial singularity. The asymptotic behavior of scale factor is determined by H±/T±. Note that 
phases III~ and III + , which are continuously connected at r = 0, are not distinguished in 
M. Since scale factor of phases III^ 1 vanishes at r = 0, they are divided into two phases in 
this paper. 
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Each region included in figure 4 and figure 5 has two branches and each branch defines 
a different phase. For example, the earlier asymptotic behavior of VII~ phase (figure 4a) 
and the later asymptotic behavior of VII + phase (figure 4b) are determined by H±/T±, the 
later behavior of VII~ phase and the earlier behavior of VII + phase are determined by the 
behavior of scale factor near r — > ±0. The phases, which have no initial singularity, are + 
phases in figure 4. 

In figure 5, phases XI and XJ/ ± are divided by V which is given by eq.(P9"D at r = 0; 
< T < 1 for phases XJ ± and V > 1 for phases XII . 



IV. DISCUSSION AND CONCLUSION 

In this paper we studied the effect of the matter in the Brans-Dicke cosmology. This was 
motivated @ by the string cosmology || with gas of solitonic p-brane by treating them as 
a perfect fluid type matter in Brans-Dicke theory. In Brans-Dicke theory, the matter has no 
dilaton coupling. From the string theory point of view, this means that the matter has R-R 
charge. So the matter considered here corresponds to D-brane gas in string theory context. 
With this matter, we found exact cosmological solutions for any Brans-Dicke parameter uo 
and for general constant 7 and classified all possible phases of the solutions according to the 
parameters involved. There are new two phases XII different from XI^ for the behavior 
of scale factor at r = 0. So the number of total phases is 19 and some of them have no 
initial singularity. We studied all the phases of cosmology numerically and gave the figures 
for the time evolution of the scale factor. 



Recently, ref. [ I !| argued that Holographic principle in the presence of a cosmological 
constant might imply the absence of initial singularity. In the Brans-Dicke cosmology, we 
do find some solutions avoiding the initial singularity. However, when we regard the Brans- 
Dicke theory as a string cosmology, we might ask whether there are solutions which resolve 
the initial singularity and the graceful exit problem at the same time. However, the cosmo- 
logical constant in Brans-Dicke theory is not a cosmological constant in string theory where 
the cosmological constant couples to the dilaton. To discuss the problem in our frame work, 
we have to consider the matter coupled with dilaton. We will discuss the problem in later 
publications. 
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FIGURES 




FIG. 1. Phase diagram of 11 regions 
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FIG. 2. Phase diagram, 12 regions 
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(a) phase / (b) phase II (c) phase III (d) phase III + 




(e) phase IV (f) phase V (g) phase VI 



FIG. 3. The behavior of the scale from phase / to VI 
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(a) phase VIT (b) phase VII+ (c) phase VIII- (d) phase VIII+ 




(e) phase IX (/) phase IX + (g) phase X (h) phase X + 

FIG. 4. The behavior of scale factor from phase VII~ to X + 
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(a) phase XI (b) phase XI + (c) phase XII (d) phase XII + 

FIG. 5. The behavior of scale factor from phase XI~ to XII + 
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